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Abstract 

An unresolved problem in J ftp phenomenology is a systematic understanding of the differential 
photoproduction cross section, da/dz[y + p — ► J/ip + X], where z = E^/E^ in the proton rest 
frame. In the non-relativistic QCD (NRQCD) factorization formalism, fixed-order perturbative 
calculations of color-octet mechanisms suffer from large perturbative and nonperturbative correc- 
tions that grow rapidly in the endpoint region, z — ► 1. In this paper, NRQCD and soft collinear 



effective theory are combined to resum these large corrections to the color-octet photoproduction 



cross section. We derive a factorization theorem for the endpoint differential cross section involving 



the parton distribution function and the color-octet J/ip shape functions. A one loop matching 
calculation explicitly confirms our factorization theorem at next-to-leading order. Large perturba- 



tive corrections are resummed using the renormalization group. The calculation of the color-octet 
contribution to da/dz is in qualitative agreement with data. Quantitative tests of the universality 
of color-octet matrix elements require improved knowledge of shape functions entering these cal- 
culations as well as resummation of the color-singlet contribution which accounts for much of the 



total cross section and also peaks near the endpoint. 
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I. INTRODUCTION 



Our current understanding of the production of heavy quarkonia is based on Non- 
Relativistic Quantum Chromodynamics (NRQCD) 0, 0, 0|, an effective theory for bound 
states of two or more heavy quarks. The NRQCD factorization formalism solves im- 
portant theoretical and phenomenological problems in quarkonium production and decay. 
Early color-singlet model calculations of Xc decay were plagued by infrared divergences j^. 
NRQCD solves this problem by providing a generalized factorization theorem which allows 
for infrared safe calculations of inclusive production and decay rates This formalism 
incorporates nonperturbative corrections to the color-singlet model, including color-octet 
decay and production mechanisms. Color-octet production mechanisms are necessary for 
understanding the production of J/ip at large transverse momentum at the Fermilab Teva- 
tron 0,0]. However, the polarization of the observed J/ip remains poorly understood 

The color-octet contribution to J/ip photoproduction gives a large enhancement to the 
cross section near the kinematic endpoint defined by z — > 1 , where z = / in the proton 
rest frame 0, Ref. 0] proposed that this peak constitutes a distinct signal for a color- 



octet contribution to photoproduction, and observed that such a signal is in contradiction 
with the experimental data which do not exhibit a peak. A crucial test of the NRQCD 
factorization theorem is verifying the universality of the color-octet matrix elements which 
enter into a variety of different J/ip production processes, and the unobserved excess of J/ip 
at the photoproduction endpoint could be interpreted as a failure of universality. 

Color-octet matrix elements were first fit to hadropro duction data from the Tevatron 
and have been refitted a number of times 11, 0, 0, w]| with general agreement among 



the various results. In general the extraction of the color-octet matrix elements has fairly 
large theoretical errors associated with it. The errors are particularly large on the extracted 
values of the color-octet matrix elements (Og ^Sq)) and (Og^( 3 Pj)), which dominate the 
color-octet contribution to photoproduction, where a variation of almost an order of mag- 
nitude is typical. These large errors are thought to be due to higher order QCD corrections 
and as a result there has been an attempt to incorporate a subset of higher order contri- 
butions through the use of Monte Carlo methods and a factorization approach [16]. 
These approaches extract values for color-octet matrix elements up to 5 times smaller than 
the fixed-order calculations. Unfortunately, the rise in the fixed-order calculation of the 
color-octet photoproduction differential cross section is so dramatic that even the large un- 
certainties in the color-octet matrix elements cannot account for the discrepancy. 

However, as one approaches the endpoint region, z — > 1, the fixed-order calculation of 
Refs. 0, [H| is invalidated by large perturbative and nonperturbative corrections. Specifi- 
cally, perturbative corrections of the form 

Jn m (l-z) 

a" i , m < 2n - 1 

1- z ~ 

are the source of the growth of the cross section as z — > 1 jl?], E3 ■ These corrections must 
be resummed to all orders to make sensible comparison with data. There are also nonper- 
turbative corrections scaling like v 2n /(l — z) n , where v is the typical velocity of the cc within 
the J/ip |l9|. These corrections invalidate the NRQCD expansion near the kinematic end- 
point. Ref. [20] addressed the latter issue by recalculating J/ip photoproduction including a 
nonperturbative shape function which resums these corrections. The shape function tames 
perturbative endpoint divergences and results in a differential cross section that peaks near 
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z ~ 0.9. However, the peak in the differential cross section is too narrow to be compati- 
ble with data, and the authors conclude that a resummation of the singular perturbative 
contributions must also be carried out. 



Similar issues arise when analyzing the production of J/ip in e + e collisions. In Ref. [21 



a color-octet factorization theorem for e + + e — > J /ip + X was derived using a combination 
of NRQCD and soft collinear effective theory (SCET) H @. This approach 
resums nonperturbative and perturbative corrections that are enhanced near the kinematic 
endpoint. An important conclusion of Ref. [2JJ is that both the resummed perturbative 
corrections and the resummed non-perturbative corrections are necessary because the two 
effects act constructively to significantly broaden the spectrum. 

In this paper we apply NRQCD and SCET to photoproduction at the endpoint of the 
spectrum. We derive a new factorization theorem for color-octet J j ip photoproduction which 
resums endpoint corrections for this process. The factorization theorem in this case is more 
complicated than that for e + + — ► J/ip + X because there is a parton in the initial 
state hadron. To derive the factorization theorem, we first match QCD onto SCETj (in 
which collinear partons have off-shellness 0( y/ MAqcb)) at the scale M = 2m c , and then 
match SCETj onto SCETn (where collinear partons have off-shellness 0(Aq CD )) |26[ at an 
intermediate scale, My/1 — z. At each stage we check that the effective field theory correctly 
reproduces the infrared physics of the previous theory to ensure that large logarithms are 
correctly resummed. Similarly, evolution of the renormalization group equations (RGE) is 
carried out in two stages. The first stage of evolution, from the scale M to the intermediate 
scale, My/1 — z, is performed in SCETi , while the second stage of evolution, from My/1 — z 
to M(l — z), involves SCETn running. 

The factorization theorem is novel in that the J/ip are not required to be produced 
with large transverse momentum, p±, with respect to the photon-proton beam axis. Most 
perturbative analyses of J/ip photoproduction assume that p± of 0(1 GeV) or larger is 
required for pQCD to be applicable to this process. For this reason, most phenomenological 
and experimental analyses impose a p± cut on the data. Such a cut is inapproriate for 
application of the factorization theorem and the resummed cross sections derived in this 
paper. 

Our final result for the color-octet contribution to da/dz exhibits a spectrum that is 
significantly broadened and has a peak reduced in height relative to including only a shape 
function or only resumming the perturbative corrections, as anticipated. The inclusion of 
both nonperturbative and perturbative corrections significantly broadens the spectrum and 
is in qualitative agreement with data. However, fixed order color-singlet calculations are 
consistent with existing data on photoproduction so if these calculations are naively com- 
bined with our calculation, the color-octet matrix elements in photoproduction need to be 
roughly an order of magnitude smaller than those extracted from fixed-order analyses of 
Tevatron data. However, there are several sources of uncertainty which make it prema- 
ture to attempt quantitative extraction of color-octet matrix elements using our calculation. 
Currently, the shape functions appearing in the photoproduction calculation have not been 
precisely determined. Also, consistency requires a similar resummation of endpoint effects 
in color-singlet photoproduction which accounts for much of the total photoproduction cross 
section. Finally, all of the existing data has cuts on p± and/or cuts on diffractive contri- 
butions which are inappropriate for our analysis. Since the resummation is expected to 
suppress the color-singlet contribution, and the cuts tends to reduce the experimental cross 
section at largest values of z, a future analysis which takes these two effects into account 
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could find more room for a color-octet contribution. 

This paper is organized as follows. In Section II, we derive the factorization theorem for 
the color-octet J /if) photoproduction cross section near the endpoint, using the two-stage 
matching procedure discussed earlier. In Section III, we obtain analytic expressions for the 
leading logarithmically enhanced corrections at next-to-leading order (NLO) in QCD. In 
Section IV, we perform the matching calculations. This section is broken up into two sub- 
sections. In the first subsection we match the SCETi current onto the QCD amplitude for 
7 + g — > cc at one loop and show that large logarithms are minimized at the scale M. In the 
second subsection we first calculate the NLO cross section in SCETi . Correct evaluation of 
the SCETi cross section requires care in avoiding the double counting of usoft modes and 
collinear modes with vanishing label momentum, the so-called "zero-bin" modes, as recently 
discussed in Ref. |27[. Taking this subtlety into account is necessary to demonstrate that 
SCETj reproduces the large endpoint corrections of the NLO QCD calculation extracted 
in section III. We then match the SCETi result onto SCETn, which determines the inter- 
mediate matching scale. In Section V, we perform the RGE evolution which resums large 
perturbative corrections. In Section VI, we discuss the impact of the resummed color-octet 
differential cross section on the phenomenology of J / ip photoproduction. A brief summary 
is provided in Section VII. Appendix A derives formulae that are useful for the extracting 
the large NLO corrections in Section III, and Appendix B provides some details about the 
zero-bin subtractions. Preliminary results describing some of the work here were presented 
at the Ringberg Workshop on New Trends in HERA Physics 2005 |28| . 

II. FACTORIZATION 

In this section we derive a factorization theorem valid for J/ifj photoproduction near 
the kinematic endpoint, which is defined by z ~ 1, where z = p^-pp/p-y-pp, and p^, pp 
and p 1 are the four-momenta of the J/ip, the initial state proton, and initial state photon, 
respectively. In the proton rest frame z = E^/E^. We begin by showing that in the limit 
z — > 1, NRQCD factorization breaks down and SCET is required to obtain the appropriate 
factorization formula. 

In the proton-photon center-of-mass frame 



where s = (p 7 + pp) 2 , = (1, 0, 0, —1), = (1, 0, 0, 1), M — 2m c , and v M and fc M are the 
4- velocity of the J/ip and the residual momentum of the cc pair in the J/if>, respectively. In 
terms of the scaling variable, z, the J/ip velocity is 




p P = 



(1) 



z*fs 



"+p'i + 



2 

mi 



(2) 



2 



n 



2z^/s 



■rv 



where m 2 , 




M 



■pI - p . 



(3) 
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In NRQCD k^ is typically dropped at lowest order, and terms with powers of k^ are matched 
onto NRQCD operators with derivatives that are higher order in the v expansion. When the 
kinematics of a quarkonium production or decay process become sensitive to any components 
of k^ such that this expansion breaks down, resummation of NRQCD operators into a shape 
function is required Il9|. 

In the rest frame of the J/ip all components of k^ scale as Aqcd- in the proton-photon 
center-of-mass frame, the components of k^ scale as 

W = (n ■ k, n ■ k, k ± ) ~ Aqcd l) , (4) 

where we have neglected 0(p±/ a/s) corrections which are small in the endpoint region. The 
n ■ k component of k^ is enhanced by yfs/M, while the n • k is suppressed by the same 
amount. Computing p\ and keeping only the n ■ k component of k^ we find 

, / M \ ( Mm 2 \ M 2 „ ^( Mm 2 \ , 
p 2 = s 1 - — z i -± - _ p 2 -yfeil -A )n ■ k + ... . 5 



In the endpoint region, the transverse momentum of the J/ip pair is of order v/I(l — z) ~ 
a/Aqcd^ ~ 1 GeV, where is the partonic center-of-mass energy which is of order M 
in the endpoint region. Therefore, the term in Eq. (JHJ proportional to is unimportant. 
Furthermore, for HERA kinematics, yfs ~ 100 GeV, so Mm\/ \szM^) ~ M 2 /s ~ 10~ 3 for 
z ~ 1. The last term in Eq. (JSJ) is as important as the first term when n-k ~ y/s(l—Mz/M^), 
which occurs when z ~ {M^ — Aqcd)/^ ~ 1- in this region NRQCD factorization breaks 
down. ^From the point of view of SCET, a new factorization theorem is required when the 
final state particles recoiling against the J/tp are jet-like and have to be described by SCET 
collinear fields rather than integrated out as in conventional NRQCD factorization. The 
final state is jet-like when n ■ px ^> \fpx-, which is equivalent to (1 — Mz/M^) <C 1. This 
leads to the conclusion that as z — > 1, SCET is required for J /if) photoproduction. 

We derive the factorization formula for the photoproduction cross section in the endpoint 
region in two steps. First we match the QCD amplitude for 7 + g — > cc onto the SCETi 
current: 

J"(x) = J d 3 u; e mv ' u) - x C^)J a (uj, x) , (6) 

where 

d 3 u; = J du)d 2 tu±, (7) 

with uj^ = um fl /2 + The leading order contribution from the color-octet 1 5'o current is 

J a (u,x) = [^{V-^Bfx-plix), (8) 
where B" is the gauge invariant collinear gluon field defined by 

B1 = ±W^n+9s(Al q )±)W, 
9s 

and if) p and Xp are the NRQCD fields for the heavy quarks and antiquarks, respectively. 
The color-octet 3 Pj current is 

J^(w, x) = A ■ p^A . - w)Bix- P ] , (9) 
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where p a = A^p* /M, A 11 " is the boost matrix from the lab frame to the cc rest frame, and 
-pM _ p n M 1 2 + is the operator which projects out label momentum. We have defined 

6®{V -u) = 8{V - oo)5 {2) {V ± - u ± ) . (10) 



for the 1 S^ S) and 3 Pj 8) channels |2jJ 



Tree level matching of QCD onto SCETi determines the leading order Wilson coefficients 

Q 

{ } ~ M C± ' 
c ^ ( 3p(8) ) = -i4eeMM) ^ f + ^ _ (n) 



The second step is to match the differential cross section from SCETi onto SCETn. For 

.(a" 



clarity we will treat the contribution explicitly and state the final result for the 3 P < - 8/! 



channel. The SCETi differential cross section is 
da f f — Cl.,C£ 



X n ,X S 

\ Y^M^i^ii °)\ J /^ + X s + X n )(J/4, + X S + X n \J a (oj 2 , 0)\p P ). (12) 



X 

2 

spin 



where p n is the total momentum of the final state collinear particles denoted by X n , and 
p s is the total momentum of the final state ultra-soft (usoft) particles, X s . The momentum 
conserving delta function can be separated into light-cone coordinates: 



<$ (4) (P 7 +Pp ~Pi> ~Ps -Pn) = 2 5(Vs(l - z) - n -p n - n ■ p s ) sfy/s - 



n-Pn 



X 



5 {2) (P±+Pn±), (13) 



where p± is the J/ip transverse momentum. The usoft momentum is a subleading contri- 
bution in the last two terms on the right- hand- side and has been dropped. As pointed out 
earlier, in the endpoint region \p±\ cMso from here on we let m\ — > M|. 

We factor the collinear degrees of freedom from the soft through a field redefinition which 
decouples usoft and collinear fields in the SCETi Lagrangian 

Bl(x) -> Y(x)Bf a (x)Y\x). (14) 

Here Y(x) is a path ordered exponential of usoft gluon fields extending from — oo to x. The 
above field redefinition also shifts the collinear fields in the out state in such a way that 
Y(x) — > Y(x) where Y(x) is a path ordered exponential extending from x to oo [29^ 

Y(x) = P exp ^ig J ds n • A us (sn + x)\ . (15) 

As a consequence the matrix elements in Eq. (J12j) can be factored into separate collinear 
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and usoft pieces. The cross section in its factored form is 



x 2 J2(J>pMT A B1(0)} 5 (3) (^ - ^ 2 )|X n )(X„|5( 3 )(P - Wl ) Tr^E^O)] |p P > 

spin 

x (01x^^^(0)17/^ + X s )(J/^ + X s |^y T B F t x _ p (0)|0), (16) 

where we have used conservation of the label momentum to replace the final state collinear 
momenta with uo and have dropped the (0) superscript on the collinear fields. The momen- 
tum components n-p n ~ 0(A QCD ) are of the same size as the usoft momentum components. 
This remaining delta-function can be expressed as an integral over an exponential 

/dec r % 
-^-exp|^-(v/i(l - z)-n-p n - n -p s )x~J , (17) 

which can be pulled into the factored matrix elements where it shifts fields from the origin 
to the point x~ . The explicit dependence on p n and p s disappears and the sums over 
final collinear and usoft states can be performed using completeness. Using d z p^j^2E^) = 
dzd 2 p±/(2z) gives 

dzd 2 p ± lQsz V M ) J + J 2 V 2 z^Ts) \ P± 2 J 

x I ^MTr^Or-) 5^(V + - ui + )S ±1 ,(0)] \p P ) 

spin 

x (0| X L p ^T A yt^ p , (a; -)^ty T A : pt x _ p(0 )| ) ) (18) 

where we simplified the collinear matrix element by projecting onto a collinear operator in 
a color-singlet configuration. The expression above is further simplified by integrating over 
LO+ and p± 

da —1 ( 2ee c g s \ f dx~ 



dz 8sz \ M 



e^( 1 -^"(o| x L p ,yr A yt^ p/(x - )w ty T A^t x _ p(0) | 0) 

1^ r / _ 2Ml \ 

x 2 E<Pp|1t[SI(x-) 5[V - -±\BU0)] \pp) , 

spin ^ V / 

~ n l 2eec9s \ (O^CSo))M J dk + S^ 1So \-V~s(l -z) + k + )J P {k + ) . (19) 



2sz \ M 



The cross section is expressed as a convolution of a shape function, S ( - 8,ls °\ and a jet function, 
Jp, that are defined as follows: 

7 4tt 2M(Of (^o)) 



spin V / 

We have set z — > 1 inside the matrix element in the definition of the jet function since the 
matrix element is smooth as z — > 1 and leads to no large logs. We will also set z — > 1 in 
the prefactor of the cross section appearing in Eq. (|19p. The shape function is normalized 
so that / dk+S^ 1Su \k+) = 1. 

Taking moments of the cross section with respect to z, = dzz N da/dz, and con- 
sidering the large N limit, we will see below that is the product of moments of the 
shape function and jet function. Large logs of M/N are contained in moments of the shape 
function and large logs of Mj \[N are contained in the moments of J7p, so the two scales are 
separated in Eq. ()19|). However, the jet function still depends two scales, the perturbative 
scale, Mj y/~N, as well as the long distance scale, Aqcd- Dependence on the latter scale arises 
because the matrix element is taken between proton states. In this respect, the jet function 
that appears in Eq. (|T9^) is quite different from the perturbatively calculable jet function 



that appears in the endpoint resummation of e + + e _ — > J ftp + X [21 . 

Since the jet function, J7p, contains both perturbative and nonperturbative scales, we can 
perform a further factorization on this matrix element. Because the external gluons in the 
proton have off-shellness O(Aqcd), the matrix element should be evaluated in SCETn rather 
than SCETi. The factorization can be thought of as arising from matching the nonlocal 
SCETj operator in Eq. ([2*T|) onto local operators in SCETn. Intuitively we expect the 
result to be a convolution of short distance coefficient which is perturbatively calculable and 
contains large logarithms of Mj y/~N with a parton distribution function for the gluon in the 
proton. 

To demonstrate this we define the variable, p = M^/s, and take moments of Jp with 
respect to p: 



1 d PP N Mk+) = fd PP N [ d ^ k+ y- 



^0 



X^(pp|Tr[^(y-) S(P - 2p^)B ±u (0)} \p P ) 

spin 



An 2 

spin 



J2(Pp\^[Bl(y-) {2V - )N+1 B ±u (0)]\Pp) ■ (22) 



Now we apply the operator product expansion (OPE) to the matrix element, which is justi- 
fied when matching onto SCETn because the momentum k + ~ 0(M/ VN) is parametrically 
larger than A QC d: 

fdpp N J P {k+) = I ^e^Wr)^^!^^-^^^]^) 

' spin 

= -±C(N,k + ) fd^fMO, (23) 
4S Jo 



where 



C(N, k + ) = J- S I ^ e ^ k+ y~C(N, y~) , (24) 



S 



is dimensionless. Therefore the moments of the jet function are proportional to the moments 
of the structure function for the gluon in the proton, whose definition in SCET is [30] 

ij>>l [Tr{BU0)6(V + -u, + )B ±1/ (0)}]\p P ) = _^/, /p (_^) . (25) 

spin 

The result in Eq. (l2~3|l is easily seen to be equivalent to the following convolution: 

Jp(k + ) = ^ / f °n (f > , (26) 

where the function Cu(p/£, k + ) can be obtained from the coefficients in the OPE, C(N, k + ), 
by inverse Mellin transform. The resulting expression for the cross section is 

Y z = aoP I dk + S^\-V~s(l-z) + k + ) ffc n (£,k+)f g/P {t), (27) 



where 



a„ = ^(0 8 J '*(%)>, 

and we have set p = Am? c / s. The factorized form for the color-octet P-wave contribution to 
the cross section is obtained by making the replacement 

(O^CSo)) - ^<0 8 JA/ fPo)> , (28) 

and replacing S^ 8 ' ls °^ with S^ 8 ' 3p °\ The shape functions in each channel are normalized the 
same way. Eq. (|2*7j) is modified once higher order corrections to the matching coefficients in 
Eq. (jllj) are included. We parametrize these corrections as follows: 

C^fS®) = C v {p) ~ 2eec ^ {M) eT (29) 

and similarly for 3 Pj 8 ^ channels. Including these corrections gives the final form of the 
factorization theorem: 

da 

dz 



a p\C v (p)\ 2 jdk + S^ s "\-^(l-z) + k + ) J^Cu^k+^fg/pit). (30) 



This factorization theorem for endpoint photoproduction of J/ip is the main result of this 
paper. For the remainder of this section we set Cy(p) = 1. In section V, where large 
logarithms are resummed, an important step is evolving Cy(p) from the scale p = M, 
where QCD is matched onto SCETi, to the scale p = M/y/N, where SCETj is matched 
onto SCETn. 

Our result should reproduce previous results for photoproduction in the appropriate lim- 
its. To lowest order in a s the coefficients in the OPE are C(N,y~) = 1, which is equivalent 
to 

C„(£,kA=Vt6{k + )s(l-£\ . (31) 
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Inserting this into Eq. (|27|) yields 



Tz = ^^(^(^0)) VSS^(-VS(l - z)) MP) ■ (32) 

To lowest order in v 2 , \/s S^ 8, So ^(—\/s(l — z)) — > 5(1 — z) which is easily seen to reproduce 
the tree level calculation of Ref . [3l| . 

Finally, it is useful to derive an expression for the moments of the cross section. Defining 
k + = t/s(u — z), the normalized differential cross section is 

ig =/ ,jr , *,s*-« ( „ )J r l | c „(£, B _ I ) /i/ , ({) , (33) 

where S^' ls °^{u) = y / i5 ,( - 8 ' l5o - ) (— y/s(l — u)) and we have rescaled second argument of Cu- 
lt is now straightforward to take moments of the cross section and show that, for large N, 
the normalized moments of the cross section factorize: 



where 




We will show in Sec. IIVI that the moments of the QCD cross section at lowest nontrivial 
order factorize in a manner consistent with Eq. ([34)1. providing additional evidence for the 
factorization theorem. 



III. EXTRACTING LARGE LOGS AT NLO 

A calculation of the leading color-octet contribution to the forward cross section for J/ ip 
photoproduction appeared in Ref. j3l|. The NLO calculation of the color-octet contribution 
to da/dz,z 7^ 1, was carried out in Refs. but these papers presented only numerical 

results. Refs. [HI 12] presented analytic results for the total partonic cross section. In this 
section we obtain analytic expressions for the leading logarithmically enhanced corrections 
to the color-octet differential cross section in NRQCD. By logarithmically enhanced we mean 
corrections of the form: 

1 ^ + ^ J l* + - Y ' oc a, 0^1) , a, (-U . (36) 



cr dz \ 1 — z 

When moments of the cross section are taken, dzz N da/dz, the distributions in Eq. (|3fi)l 
give rise to terms oc a s ln 2 iV and a s lniV, respectively. These are the logarithmically enhanced 
terms we wish to resum to all orders in perturbation theory. Terms in the differential cross 
section less singular in 1 — z give contributions to the moments which are suppressed by 
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powers of N, e.g. J^dzz N \n(l - z) = kiN/N + 0(1/N 2 ), where N = Ne~< E , and will be 
ignored. 

The analytic result for the logarithmically enhanced corrections is helpful for checking 
our factorization theorem for the resummed cross section. Furthermore, in our derivation of 
the factorization theorem, QCD currents are first matched onto SCETi operators, then the 
cross section in SCETi is matched onto SCETn. In both matching calculations, we wish to 
verify that no large logarithms appear in the matching coefficients so we are sure that large 
logarithms are properly resummed using the RGEs. The condition that the large logs cancel 
in the matching from SCETi and SCETn determines the scale at which SCETi is matched 
onto SCET n . 

The hadronic cross section is obtained by convoluting the partonic cross section with a 
parton distribution function (pdf): 

dz ^-^ l n x 2 \xJ dz 

Here p = M 2 / s and x = M 2 / s, where s is the center of mass energy squared in the partonic 
collision. There is a sum over parton species, i, and fym is the pdf for the proton. For 
the rest of the paper we will restrict ourselves to the dominant production process which is 
initiated by a gluon. In the HERA collider J/if> photoproduction experiments ^/s ~ 100 GeV 
so p ~ 10 -3 . Though we cannot set p to zero inside the argument of the pdf, we can set 
p = in the limit of the x integration since the logarithmically enhanced terms come from 
the opposite end of the integral. Calculation of the partonic cross section is simplified by 
making this approximation, so we will use the standard distributional identity 



'l-xr^ = --5(l-x)+l — ) - e (-5— +0(e 2 ), (38) 



instead of ^ 

(1 - x)- 1 " 6 = -^6(1 - x) + (-^—) - e ~ X) ) +0(e 2 ), (39) 



x J \ 1 — X 



where (3 = y/1 — p and the p-distributions are defined by 

dx[t(x)] p f(x)= I dxt{x)[f{x)-f{\)\. (40) 



By using Eq. rather than Eq. ()39|) we omit corrections proportional to log/3 ~ 10~ 3 
which can be safely neglected. 

The leading order contribution to the production of J/ip via color-octet mechanisms is 
the two-to-one process, 'jg — > cc^( 2S+1 Lj), where one of the two Feynman diagrams is 
depicted in Fig. [TJ This gives a contribution to the partonic cross section [3l| 



da LO 7f 3 a s ael 
dz Arw> 



5(1 -x) 6(1 -z) 



= a 6(1 - x) 8(1 - z) . (41) 
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A 



FIG. 1: One of two leading order QCD diagram for 7 g — > cc^ ( 2S+1 Lj) 
is obtained by reversing the direction of the fermion line 



The other Feynman graph 




FIG. 2: An example of a NLO QCD diagram contributing to photoproduction. 



An example of a next-to-leading order real emission diagram that contributes to the pho- 
toproduction cross section for z 7^ 1 is shown in Fig. |2j It is these graphs that contribute 
the large logarithms at NLO we wish to resum. Infrared (IR) divergences in these graphs 
will be regulated using dimensional regularization and the IR divergences are cancelled by 
virtual corrections to the LO graphs, Fig. ^ or by the parton distribution function. 

The partonic cross section obtained from summing all the NLO real emission diagrams 
in the color-octet 1 So channel is 

d<r real C A a 8 (A-Kn 2 YxHl-z)-Hz-x)- € ... . 

Sr J, , M(x,z) , (42) 



dz n \ M 2 J T[l 
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where 



M(x,z) 



x 2 (z — x)[(l — z + z 2 ) 2 + x(x — z)(x(x — z) + 2z 2 )] 



1 



X) 



[1 + x 



z)z 2 



(43) 



The phase space is in D = 4 — 2e dimensions. The logarithmically enhanced terms come 
from the limit z — > 1. Expanding M(x, z) about z = 1 we find 



M(x,z) 



X + X 



2\2 



[l-x)(l-z] 



+ 



-2 + 5x- 6x 2 + 5x 3 - 6x 4 + 5x 5 - 2x 6 
x(l — x) 2 



+ 0(1 - z) . (44) 



Terms 0(1 — z) and higher do not give logarithmically enhanced contributions, even after 
performing the x integration, and therefore can be dropped. Furthermore, 



dx- 



(X 



X) 



dx- 



1-x) 



oc ln(l — z) 



(45) 



so terms proportional to 1/(1 



give logarithimically enhanced corrections while terms 



proportional to 1/(1 — x) can be dropped. Thus the second term on the right hand side of 
Eq. (J44j) can be expanded around x — 1, 



da 



real 



dz 



0- o - 



7T 



M 2 



xHl- z) 



z — X) 



~M 2 



r[l - e] 

x e (l — z)~ e (z — x) 



2\2 



[1 — X + X 

[(l- x )(l-z) 
2C A 



1 - z 



l-x) 



(l-x) 2 

+ . 



+ ... 



(46) 



where the ellipsis represents terms that do not give logarithmically enhanced corrections and 
P gg {x) is the related to the real emission part of the gluon splitting function: 



2C A 



x 



+ 



X 



+ x(l 



x 



X 



X 



The gluon splitting function, V gg (x), will appear in what follows and is defined by: 



^gg( x ) = P 9 g{x) + b 5{l-x) 



Pggi x ) 



bo 



2C A 
11 



x 



[l-x) 
2. 



+ 



x 



+ x(l — x) 



X 



6 



■C A - -T F n f 



(47) 



In Eq. ()46|) the first term in square brackets contains the collinear divergence, as is easily 
seen by noting that 1 — z = (1 — x)(l — cos 0)/2, where cos 9 is the angle between initial and 
final state gluon momentum in the parton center of mass frame. 

Though we have derived Eq. (|46|) for the 1 Sq channel, the result is in fact universal and 
holds for the color-octet channels 1 S'o, 3 -Po,2 regardless of the angular momentum quantum 
numbers of the final state cc. QCD factorization theorems guarantee that the term with 
the 1/(1 — z) pole in the real emission cross section has the same form for all color-octet 
production channels. The universality of the double pole, 1/(1— a;) 2 , is verified by calculating 
QCD diagrams in the soft gluon approximation, which is valid in the limit x — > 1 (32| . 
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To obtain the logarithmically enhanced corrections we need to perform the integration 
over x. At first sight, deriving an analytic expression appears to be difficult since the pdf 
contains x dependence which is not known in closed form. However, it is possible to extract 
the leading corrections analytically using the following distributional identities, derived in 
Appendix A: 



:i-z) 



-l-e 



dx 



Z - X 



X 



5(1 -z) 
1 

+ " 



1 - z 



1 

2e 2 12 
i 

dx 



9(1) 

l 



dx 



x 



g(x) + 



dx 



g(x) - g{ 



dx- 



[z — x) 



:i-z) 



X) 



-^(i 



V i - z 
1 



1 - z 



9(1) 



ln(l — x) 



x 



(48) 
g(x 



(49) 



where g(x) is arbitrary. On the right hand side we only keep terms proportional to singular 
distributions in 1 — z. This includes the singular distributions in Eq. which contribute 
to the large logarithms in moment space as well as terms proportional to 5(1 — z). These 
terms contain IR divergences which cancel against virtual graphs or are absorbed into the 
pdf. 

The NLO result for the part of the partonic real emission contribution to da/dz that is 
singular as z — > 1 is 



da 



real 



dz 



C A a s 
a x 



7T 



(50) 



5(1- z) {5(1- x) 



+ (1 — X + X 

1 



2\2 



2e 2 
1 

e 



2e 2e 



P 2 



M 2 



P 2 



In 



— 2 

H x 



x 



+ 



1 - z 



-5(1 -x) + (1-x + x 2 ) 2 



1 — X 



M 2 



ln(l — x) 



P 2 



M 2 



7[- 



-5(l-x) 



/ ln(l - z) 



where we define /i 2 = 47r/i 2 e lE . If we take moments of this cross section with respect to 
the variable z we find 



= 5(1 - x) x 



7T 

1 11, 

2T 2+ Ye + Ye ln 



o:. s 



(51) 



M 2 



-In 



+ lln 2 

4 



fx 2 xN 
M 2 



M 2 

+ ^-X(l - X)Pgg(X 



1, 

+ -In 

2 



M 2 



+ ln(A)-lln 2 (A)-^ 



/ln(l 



x) 



2tt 



V i 



X 



This expression for the partonic QCD cross section is needed for the matching calculations 
in the next section. 
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IV. MATCHING 



In this section, we match QCD onto SCETi at the scale ll = M and then SCETi onto 
SCETn at the scale Jx = Mj Vn. At the first stage we must verify that SCETi reproduces 
the IR behavior of QCD, i.e. the SCETi real emission diagrams must reproduce the terms 
in Eq. (|4*U|) which dominate as z — > 1. At the second stage of matching, the SCETi cross sec- 
tion is matched on the SCETn cross section, which takes on the factorized form of Eq. (|27j). 
These calculations are simplest if dimensional regularization is used to regulate both ul- 
traviolet (UV) and IR divergences. The matching calculations are important for verifying 
the factorization theorem and determining the QCD-SCETi and SCETi -SCETn matching 
scales, which give the boundary conditions for renormalization group evolution. 



A. Matching QCD onto SCETi 

In the first step we match the QCD amplitude onto the SCETi current at one loop. The 
one-loop QCD result for the virtual correction can be found in Ref. |32T |. 



^ = ((to + a^)S(l - x)5(l 



(52) 



a 



(V) 



a 



7T~ 



Cf — rCi Ca -tt 



17 



6ei R / 3ei R 



n f T F + 2D 



o 



where Cf = 4/3, Ca = 3, Tp = 1/2 are SU(3) group theory factors, nj is the number of 
light flavors, v is the relative velocity of the c and c in the rest frame of the cc pair, and 



47T// 



M 2 



So 



D 3 



Po 



C F 
C f 



r(i + e) 

+ y )+C A [^ + - 

7 7T 2 

'3 + T 



+ C A 



[8] 
Pi 



-4C F + C A ( - + 



log 2 



7T 

+ T 



7T 
12 

12 

2 



(53) 



The UV divergence in Eq. (jo^jl is removed by QCD coupling constant renormalization. 
When computing the matching coefficient to this order, one must include one-loop NRQCD 
virtual corrections to the matrix element of the NRQCD color-octet production operators. 
The NRQCD corrections reproduce the Coulomb correction oc tt 2 /v in Eq. (|52jl and therefore 
this term does not appear in the matching coefficient 33]. 

In SCETi loop integrals are scaleless in dimensional regularization and therefore are zero. 
We must also subtract the contribution from the EFT counterterm which was calculated in 
Ref. E 



Z -l 



4n 



i i 



log 



M 2 



+ 



17 
67 



3e 



n f T F 



(54) 
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FIG. 3: Feynman graphs for the real emission of a collinear gluon. 



Thus the SCETi one loop result for the virtual cross section is 



do 

dz 



a \C v (p)\ 2 5(1 - z) 1 + 



2tt 



1 1 , 



\M 2 



+ 



17 

6e 



+ ¥ n f T F 



• (55) 



Taking the difference between this and Eq. (|52|) (after subtracting the UV divergence and 
dropping the Coulomb correction) we obtain the one loop matching coefficient 



C v ip) = 1 + 



47T 



/ - , . T . 2 y" 2 17 1 A* 2 71-2 



o 



(56) 



The logarithms are minimized for /i = M, which fixes the QCD-SCETj matching scale. 



B. The NLO SCETi differential cross section and matching onto SCETn 



In the second step we match the SCETi differential cross section onto the 
SCETn differential cross section which has the factored form given in Eq. (|27|). The calcu- 
lation of the SCETi differential cross section at NLO also allows us to confirm that the EFT 
reproduces the parts of the NLO QCD calculation that are singular as z — > 1. This is an 
important check on the validity of the EFT. 

The SCET diagrams for the real emission come from both collinear and usoft graphs. 
These are shown in Figs. Q and (@J), respectively. The result of the usoft graphs is the 
same as making the eikonal approximation in the full QCD diagram [32j: 



d,7 SOft 

dz 



a - 



Caol s £ e (l — z ) 



X) 



TV 



r[i 



X 



X 



(l-z){l-x) [l-xf 



(57) 



Evaluation of the collinear graphs is more subtle. A naive evaluation employing the Feynman 
rules of Ref. (2l| yields 



do co1 

dz " 



CaO-s £ e (l — z) e (z — x)~ 



71 



ri 



\2\2 



X (1 + x — z + (x 
(1 - z)(z - x)(l + x - zf 



(58) 



Note that the naive collinear contribution to the cross section has pole in z — x that is not 
present in the full QCD calculation in Eq. ()43|) . Combining the cross sections in Eq. (15 7|) 
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FIG. 4: Feynman graphs for the real emission of a usoft gluon. 



and Eq. (|58|) then expanding about z = 1 as was done earlier for Ai(x, z) does not yield the 
two leading terms in Eq. (J46|) . 

The origin of this problem is double counting of modes in certain corners of phase space 



where modes can be considered either collinear or usoft 27]. Collinear particles in SCET 
have momentum scaling as Q(A 2 ,1,A). For the bulk of the phase space integration, this 
scaling is satisfied but when x, z — > 1 simultaneously it is not. To see this note that the 
momentum of the final state gluon in the real emission diagrams is 

^VB{1 - z)^ - -■ 



k» = -Vs(l - z)n^ - -Vs(x - z)n^ + k^ , 



h\ = -§(x - z)(l - z) . (59) 
For 1 — z ~ 0(X 2 ) and x — z ~ 0(1) we see that the final state gluon momentum scales as 

(n ■ k,n ■ k, k±) ~ V^A 2 , 1, A) , 



appropriate for a collinear particle. However, the phase space integral over x ranges from 
(nearly) zero to z, so there is a corner of phase space in which x — z ~ 0(X 2 ) rather than 
0(1) and therefore 



(n-k,n-k,k ± ) ~ Vs(A 2 ,A 2 ,A 2 ). 

In this regime the scaling is appropriate for a usoft particle. 

In the SCET label formalism, double counting arises from a naive evaluation of the sums 
over collinear label momenta appearing in loops and phase space integrations. The mode 
with vanishing label momentum, the so-called "zero-bin", must be excluded from these 
sums since a mode with vanishing collinear label momentum should be considered soft. 
As discussed in Appendix B, the standard trick for converting sums over label momenta 
with integrations over residual momenta to obtain ordinary loop and phase space integrals 
actually includes the zero-bin mode, and therefore certain zero-bin subtractions must be 
made to correctly evaluate collinear diagrams. 

In our problem, the zero-bin subtraction can be implemented by including additional 
diagrams that are identical to the collinear diagrams except that the momentum of the final 
state gluon is treated as usoft and the appropriate approximations are made. This graph 
gives 



da zb C A a s x e {l- z)- e (z - x)~ e \ -x 2 



dz ° 7T r[i 



z)(z — x) 



(60) 
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The minus sign reflects the fact that we must subtract the zero-bin contribution from the 
naive evaluation of the collinear contribution. Combining the usoft, collinear, and zero-bin 
subtraction yields 



da 

dz 



C A a s x € (l 



x) 



71 



X 



r[i - e] 

'1 + x — z + (x — z) 2 ) 2 



a - 



[1 — z){z — x) 
C A a s x e (l - 



[l + x-zf 
'1 + x 



1 + 



7T 



ri 



X 



2\2 



X 



X 



X] 



z)(l-x) (l-x) 2 _ 

777 + • 



X) 



(61) 



where the ellipsis represents terms which do not contribute logartihmically enhanced correc- 
tions. We see that once the zero-bin contribution is subtracted, the spurious pole at x = z 
is removed and the leading behavior of the QCD cross section is reproduced by SCETi, as 
expected. 

To get the SCETi differential cross section to this order we must add the virtual contribu- 
tions. Loops are scaleless and therefore vanish in dimensional regularization, but the contri- 
bution from tree level graphs with the counterterm in Eq. (p>4"j) must also be included. We will 
match the moments of the partonic cross sections rather than the partonic differential cross 
sections. Combining the SCETi real emission contribution with the SCETi counterterm 
contribution, we obtain 



O"0 



Z7T 



.'/.'/ 



-5(1 -x 



(x)- 
e 

71 



a 



— XP gg (x)\li 

Z7T 



i" 2 

Ji 2 xN 
M 2 



P 2 



M 2 



Z7T 



M 2 



+ ln(iV)-iln 2 (iV)-^ 



ln(l -x) 
1 — x 



(62) 



In SCETn, the partonic cross section takes on the factorized form of Eq. (J34)) . except now 
the initial state is a gluon rather than a nucleon: 



In Eq. ()63|) . we need the one loop result for the moments of the shape function, 



(63) 



S^\N) = 1 



a,C_ 



71 



log 2 



log 



flN 



(64) 



and the one loop expression for the gluon structure function |3 



wo=*(i-o-^a) 



(65) 



Parametrizing the matching coefficient Cn as 



x 



X 



X 



-Cn ( 7, N 1 = 5 ( 1 - j 1 (1 + ctsCff(N)) + a s Cf> ( - t , N 



.(2) 



X 



(66) 
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inserting into Eq. (|63|) and expanding to 0(a s ) we find 



a 



11 

N 



0Q 



6(1 -x) 



l + a 8 Cf}(N) 



a,Ci 



71 



log 2 



ftN 



log 



pN 

~W 



--^xV gg {x) + a s CP{x,N). 
e In 



Comparing Eq. and Eq. (JBTj) we obtain the matching coefficient Gu to 0(a s ): 



C n (x,N) = 5(1 -x) ^1 + 

(x)ln 



Us-p 
'27T " K 



C A a s 


2 




71 




( p, 2 xiSi 


: ) + 




1 M 2 




-(i- 



7T 



ln(l — x) 



(67) 



(68) 



This result shows that all large logs of N in the matching coefficient vanish at the interme- 
diate scale pL = Mj y/~N, which fixes the SCETi -SCETh matching scale. At this scale there 
remains a ln(x) in the matching coefficient which is 0(1) in our power counting. 



V. RUNNING 



In this section the renormalization group is used to resum large logs of iV in the moments 
of the cross section. An inverse Mellin transform of the resummed moments is performed 
to obtain an analytic resummed formula for the differential cross section. The evolution 
is carried out in two stages. First, we run the effective theory currents in Eq. © using 
the SCETj RGEs from the scale Jx = M, where QCD is matched onto SCETj, to the scale 



p, = M/vN, where SCETi is matched onto SCETn- Second, the shape function in SCETn 
is run down to the scale p, = M/N. 

The currents that arise in photoproduction also arise in radiative decays, and the running 
was calculated in Ref. \34 



log 



Cy(Li) 

C V (M)_ 



A7iY\ 



adj 



(3la s (M) 
Bt + 2 7l 



- 1 + log y 



PI 



1 - y + ylogy - ^ log 2 y 



0o 



4r adj 

hgy-—^- 

Po 



y - 1 - log y 



(69) 



where y = a s (/x)/a s (M) and 

rf = c A , rf = c A 



-Ca 



7i = — 



0o 
4 



0o 



--)- 


5n f 


12 J 
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2 


3 A 


— nt 
3 / 



(70) 

. At leading 



The two loop cusp anomalous dimension was first calculated in Refs. 
order C V (M) = 1. 

In the second stage, only the shape function is evolved. This is possible since the soft 
and collinear sectors in SCETn are decoupled. Denoting the scales in the collinear and soft 
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sectors as « c and fi s respectively, and making the scales explicit in Eq. (jSUj). we have 
da 

dz 



pvoCUflc) I dk + S^ S °\V^(l -z)- k + ;fi s ) 

zjfg/p{^\ fic) ■ 



(71) 



In this section, we focus on the contribution to the cross section. The expression for the 
resummed cross section is easily generalized to 3 Pj channels since the evolution equations 



are identical. In matching SCETi onto SCETn we must set u s = \i c = Mj v N in order to 
minimize logarithms in the matching coefficient Cn- However, the shape function will still 
contain large logarithms of N. The shape function must be evolved from p, s = Mj V~N to 
pL s = Mj N to minimize its logarithms. 

The running of the soft function is easily carried out in moment space where the differ- 
ential cross section is a product of the moments of the shape function and moments of the 
hard coefficient, as in Eq. 



a N = pa CUfic)S {8 ' So \N ] fi 



(72) 



The anomalous dimension of the shape function can be calculated by replacing the J/ip 
projection operator V$ = aLa^ with an on-shell charm projection operator V C c since the 
renormalization of an operator is only sensitive to short distances. The tree and one gluon 
Feynman rules are the same Feynman rules as those of the shape function that appears in 
quarkonium decay, and the one loop calculation of the shape function in Ref. j34| gives 



l(N;fl) 



71 



1 - 2 log 



flN 
~M~ 



(73) 



The RGE can be solved in moment space. Combining the running in SCETi of the currents 
in Eq. (|69J) with the evolution of the shape function, we obtain the following resummed 
expression for the moments of the J/ip photoproduction cross section: 

pa e lo ^ N)9lM+92M S^ 1So \N;M/N) 



a N 



x 



1 jC II ^,N ] M/^N)f 9/P ^ 1 M/VN) 



(74) 



where 



9i(x) 
<fe(x) 



2r adj 
gr adj 



[(1 - 2 X ) log(l - 2 X ) - 2(1 - X ) log(l - X)) , 



[-log(l-2x) + 21og(l-x)] 



log(l - 2 X ) - 21og(l - x) + - log 2 (l - 2 X ) - log 2 (l 



X) 



471 2Bi 

log(l - x) + —r~ log(l 

Po 



ATI 



adj 



00 



log n [log(l - 2x) 



2x) 
log(l - x)] 



(75) 
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with n = e 7E , x = log(JV) a s (M)/3 /47r, and = (34Ci - lOC^n/ - 6C F n f )/3. 

To obtain the resummed expression for the differential cross section, da/dz, we take the 
inverse Mellin transform of the expression in Eq. (|74p. However, this is complicated since 
the scale in the pdf depends on N. To perform the inverse Mellin transform we must first 
undo the convolution in Eq. (J74j) by taking moments with respect to p: 



dpp K 1 a N (p) 



0o 



1/' 

e log(AT) S i (x)+92(x) ( ^(8. 1 S'o) / jy . / j\A 



(76) 



X 



dp p 



K 



(p 



-r-Cn ^,N; M/VN / ff/P (£; M/VN) 



; iog(Jv) ffl (x)+9 2 (x) ^(s. 1 ^)^. ^ jy. M/VN) f g /p(K; M/VN), 



where 



C n (K, N;p) 

fg/p(K) p) 



dU K f g/P (ti;p). 



(77) 



The dependence of the pdf on the scale M/ vN can be made explicit by using the evolution 
equations for the moments of the structure function. Ignoring mixing, the running of the 
pdf in moment space is given by 



(78) 



where the explicit form of a gg (K) is not needed in what follows. The leading order solution 
of Eq. (J78j) in moment space is 



l/P [K] M/VN) = f g/P (K- M) exp 



a gg( K ) 
Wo 



log(l - X) 



= f g/ p{K; M)exp [h K ( X )).{79) 



The pdf on the right hand side of this equation is evaluated at the scale M, and is therefore 
independent of N. All N dependence has been moved into the factor hx(x)- Using this 
result in Eq.([75|) we get 

cr 



Jo g (N) gi ( x )+ g2 ( x )+h K ( x )g(8,is )( N . M j N ^ c n (K, N; M/VN) f g/P (K; M) . (80) 



Since the logarithms in Cn(x, N) are minimized at the scale M/ v N, Cu(K, N; M/ v N) = 
1 + 0(a s ), where the 0(a s ) term has no logarithmically enhanced contributions. Therefore, 
we can set Cji(K, N; Mf VN) = 1, making it possible to analytically evaluate the inverse 
Mellin transforms with respect to both N and K. Using the results of Ref. 3^ to evaluate 
the inverse Mellin transform with respect N yields 



da(z, K) 
dz 



00 



1 ^s^iz/u) 
u 

-u— { (-)( 1 - //) 



du 



e lgi(l)+92(l)+h K (l) 

T[l - 9l (l) - lg[(l 



f 9 /p(K;M), 



(81) 
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where I = —a s Po/(4n) log(l — u) and g[(l) = dg\(l)/dl. Next, we eliminate the factor Hk{1) 
using the leading order solution of Eq. (|75|) again, 



e h «Vf g /p(K; M) = f g/P (K; M) exp 

Using this result in Eq. (j81|) we arrive at 
da(z, K) f 1 du 



^log(l-0 



f g/P (K;MVT^). (82) 



S {s ' 1So) (z/u) (83) 

dz J z u 

The K dependence is now entirely contained in the moments of the pdf, so the inverse Mellin 
transform with respect to K is trivial. The fully resummed differential cross section is 

£ = P*o f-S^\z/u) (84) 
dz J z u 

( d ( e «9i (0+92(0 i \ 

x | 9(1 - u) rn - - im] f " F(n M ^ u) \) ■ 



This result is unusual because the scale of the pdf varies between My/1 — z ~ ^MAq C d 
and in the convolution integral. It is not clear how to evaluate the pdf when the scale 
is below Aqcd- Another problem arises in the resummed coefficient since the functions 
Qi(V) blow up for u sufficiently close to 1. If the resummed exponent is expressed as an 
integral over a running coupling this divergence is due to the Landau pole. This problem 
commonly arises in resummed calculations and a prescription for dealing with the Landau 
pole is required for phenomenological work. In what follows we cut off the integral at 
an upper limit, u max = 0.93, which is set by the location of the Landau pole. For this 
value of u max the lowest scale at which the pdf needs to be evaluated is 800 MeV, so the 
prescription for dealing with the Landau pole in the resummed exponent also fixes the 
problem associated with the scale of the pdf. The difference between integrating to u max 
as opposed to one is formally the same order as power suppressed corrections and therefore 
can be systematically neglected [3^|. Likewise, the difference between the prescription used 
in this paper and other possible prescriptions is of order power suppressed corrections. An 
alternative approach which explicitly avoids the Landau pole is given in Ref. jloT ]. 



VI. PHENOMENOLOGY 



Before we can investigate the phenomenological consequences of our resummed cross 
section, we must determine the shape functions. Unfortunately not much is known about 
these nonperturbative functions, although they also arise in e + + e~ — > J/ijj + X and 



19l |. A fit of the shape function to Belle 41 1 and Babar ji^ data on 



electroproduction 

e + + e~ — > J /ip + X was carried out in Ref. j2l|, and we will use the parameters determined 
from that fit in our calculation. We use these shape functions for illustrative purposes 
only, as the shape functions extracted in Ref. jij are not reliable. The Belle collaboration 
observes that the J/ip cross section at y/s = 10.6 GeV is dominated by events with an 
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additional cc pair, i.e., e + + e — > J/ip + c + c + X 4^, 44]|. This fact is poorly understood 
at the present time 45|| and we do not expect the color-octet mechanism studied in this 



paper to contribute significantly to e + + e" — > J/ip + c + c + X . A sensible extraction of the 
color-octet shape function from e + + e~ — > J /if} + X requires removal of events with extra 
cc pairs, and such data is not currently available. 

For simplicity, we assume that the 1 S < ^ > and 3 -Pj 8 ^ shape functions are the same, but this 
need not be the case. The shape function model we adopt is a modified version of a model 
used in the decay of B mesons |46j |. 



f(r + ) = ^~j{x - l^-V^- 1 ) 6(x-l), x= r T , (85) 

where a and b are adjustable parameters and A = — M. This function is related to 
the shape function only in the J /if) rest frame where /(— r + ) = S(r + ). To relate the model 
above to the shape function S(u) which appears in Eq. (J8*lj) we use boost invariance 

dk + S{-^{l - z) + k + ) -> dl + S(-M(l - z) + 1+) , (86) 

where l + = (M/^)k + is the rest frame residual momentum which is O(Aqcd)- Since 
dk+ = (y/s/M)dl + we get S(-y/s(l - z) + k+) -> (M/y/s)S(-M(l - z) + /+). Finally we 
get 

S(z-u) = y/ZS(-Vs(l- z) + k + ) 
-> M^(-M(l - z) + / + ) 

= Mf(M(l - z) - /+) = Mf(M(u - z)) , (87) 

where in the J/-0 rest frame Z + = M(l — w). 
The first three moments f{r + ) are 
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m = / <ir + /(r + ) = l, rrii = / c/r + r + f(r + ) = A(b + 1), 

m 2 = dr + (r+) 2 /(r+) = A 2 Q + (6 + l) 2 ^) . (88) 

Since A ~ O(Aqcd), any choice with a ~ 6 ~ ^(1) gives the desired scaling for the moments. 
We use parameters taken from a fit to the e + + e~ — > J /if} + X data j2l|: a = 1, 6 = 2. The 
value of the first and second moments of the shape function for this choice of parameters are 
890 MeV and (985 MeV) 2 respectively. Since m c v ~ 500 MeV the moments are consistent 
with the velocity scaling rules. 

Calculations of da/dz are shown in Fig. |3J The solid line is the final result, Eq. (|84j) . 
normalized to o"o using the shape function in Eq. (|85|). The dashed line is a plot of the shape 
function alone and the dotted line includes only the perturbative resummation. Here we 
use m c = 1.4 GeV, = 100 GeV, and the CTEQ5L pdfs 0. Fig. El shows the effects 
of both the perturbative resummation and the shape function. The results are very similar 
to other resummed calculations: the perturbative resummation alone or the shape function 
alone give a cross section that is highly peaked in the endpoint region. In contrast, the 
convoluted result is much broader with the peak shifted to lower values of z. 

It is important to point out that the calculation is not valid for z very close to one. 
Though SCET is valid in the endpoint region, the calculation assumes that the final state 



23 



FIG. 5: The color-octet contribution at the endpoint normalized to ctq. The solid line is the pertur- 
bative resummation convoluted with the shape function. The dotted line is perturbative resummation 
only, and the dashed line is no resummation, but including a shape function. 

is inclusive. However, as z gets very close to one the cross section is dominated by exclusive 
final states. This occurs when M(l — z) ~ Aqcd> which is roughly z ~ 0.9. In this exclusive 
region our analysis is not valid and a very different approach is needed. 

In Fig. IH1 we show the differential cross section. The short dashed line is the color-octet 
contribution, which has a normalization set by the linear combination of color-octet matrix 
elements in Eq. 1)41)1 . We set this combination to 3 • 10~ 3 , which is an order of magnitude 
smaller than the value determined from a fit to Tevatron data. The long dashed curve is the 
color-singlet contribution, and the solid curve is the sum of singlet and octet. The singlet 
contribution is peaked in the endpoint and it is not possible to compare to data until a 
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FIG. 6: The differential cross section in nb. The resummed color-octet contribution (short dashes), 
the leading order color-singlet contribution (long dashes), and the sum (solid line). The rise in the 
differential cross section at the endpoint due to the color-singlet contribution must be resummed 
before a meaningful comparison to data can be made. 

resummation of the singlet contribution is completed. This work is in progress [i^ j. 

The analysis of Ref. indicates that the data is well described by the color-singlet 
cross section alone so it appears that there is little room for a color-octet contribution. A 
naive fit with the leading-order color-singlet calculation and the resummed color-octet cross 
section would lead to the conclusion that color-octet matrix elements are smaller than those 
obtained from fitting to Tevatron data. However there are caveats. First, large endpoint 
effects in color-singlet production also need to be resummed. While we do not expect the 
results of this resummation to be as dramatic as in color-octet production, resummation 
should also suppress color-singlet production in the endpoint region, creating more room 
for a color-octet contribution to photoproduction. Second, most of the available data has 
cuts that remove any contribution for p± < 1 GeV. In our calculation we integrated over p± 
assuming p± < M. Thus care must be taken when comparing to experimental measurements. 
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The Zeus collaboration has carried out an analysis which includes a measurement of da/dz 
with no p± cut as well as a cut of 1 GeV and 2 GeV [49] . The data clearly show that 
a p± cut decreases the differential cross section with the greatest reduction being in the 
high z bins. Thus comparing to data with no p± cut may allow for a larger color-octet 
contribution. Furthermore it is common practice to cut out the diffractive contribution 
from the data. However, our calculation assumes that the diffractive contribution is part 
of the complete sum over final states, and an honest comparison requires data without the 
diffractive contribution subtracted. The Zeus analysis of Ref. ji^] gives the percentage by 
which the cross section decreases when the diffractive contribution is removed. Again this 
fraction becomes large as the endpoint region is approached, and including the diffractive 
contribution in the experimental measurement could lead to greater room for a color-octet 
contribution. 

VII. SUMMARY 

In this paper we studied the color-octet contribution to J/ip photoproduction near the 
kinematic endpoint. As z — > 1 the usual NRQCD factorization formalism breaks down due 
to large perturbative and nonperturbative corrections. We combined SCET with NRQCD 
to derive a factorizatioin theorem for the differential cross section, da/dz, in the endpoint 
region in terms of the parton distribution function and nonperturbative color-octet shape 
functions. Large Sudakov logarithms which appear in the endpoint are resummed using the 
RGE of the effective theory. 

Since the total photoproduction cross section is largely accounted for by the color-singlet 
term in NRQCD a resummation of this contribution must be carried out before a quantita- 
tive comparison to data can be made. However, some qualitative conclusions can be drawn. 
First, we find that perturbative resummation acts constructively with the shape function 
to significantly broaden the z distribution. This important effect improves the agreement 
between the shape of the theoretical prediction and data. Second, the normalization of our 
prediction is set by a linear combination of NRQCD matrix elements which also appear in 
calculations of J/ ip production at the Tevatron. The data clearly prefer values for the color- 
octet matrix elements that are roughly an order of magnitude smaller than the central values 
extracted from Tevatron data. However, serious comparison with data requires resummation 
of large endpoint corrections to the color-singlet production cross section. Also, the assump- 
tions underlying our calculations require that we apply our results to data without cuts on 
p± and diffractive contributions that are commonly used in existing experimental analyses. 
Resumming color-singlet contribution and relaxing experimental cuts could both lead to 
more room for color-octet contributions in a complete analysis of J/ ip photoproduction near 
z->l. 
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APPENDIX A 



Here we derive Eqs. ()48|) and (jUj) . We aim to extract terms that are singular as z — > 1, 
so we will systematically drop all contributions which are regular in this limit. We start 
with Eq. fllHJ): 



'l-z) 



-l-e 



dx g(x) ■ 

1 — x 

(l-z)" 1 - 



(Al) 



g-(l) / 



^ — X) 
1 — X 



+ dx 



v z — X) 
1 — X 



-(g(x)-g(l)) 



We assume g(x) can be expanded about x — 1. To evaluate the second term on the right 
hand side of Eq. (jAl|) . we note the distributional identity in Eq. (JHHj) implies 



'l-z) 



-l-e 



/(*) 



-l-e 



(A2) 



where the ellipsis denotes terms that are not nonsingular as z — > 1 for any function f(z) 
which can be expanded about z — 1. The integral in the second term can therefore be 
evaluated with z = 1: 

(z — 



-l-e 



1 — X 



(g(x)-g(l)) 



(A3) 



l-z) 



'l-z) 



-l-e 



-l-e 



dx(l-x)- 1 - e (g(x)-g(l)) 



dx 



1 — x 



ln(l — x) 

1 — x 



We can use the distributional identity in Eq. (J3~%j) and expand in e to obtain 



'l-z 



,-l-e 



dx { \ X) ~ € (g(x)-g(l)) 
1 — x 

8(1 -z) 
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1 — x 



g(x) + I dx 



ln(l — x) 
I — x 
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l-z 



dx 



1 — x 



To evaluate the first integral on the right hand side of Eq. (|A1J) we remember that 



(l-z) 



-l-e 



dx 



(z — x) 
1 — x 



is really a distribution so we consider the double integral: 



dzf(z){l-z)- 



-l-e 



dx 



(z — x) 
1 — X 



where /(z) is a smooth function. Writing / as 



dz\f(z) -f {!)]{!- z)- x ~* 
•i 



dx 



1 — X 



+/(1) / rfz(l-z)- 1 ^ 
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1 — x 



we see that the integrand of the first term is finite as z — > 1 so we can set e — > in evaluating 
this term. The second term is straightforward to evaluate after interchanging the order of 
the x and z integrations. The result is 







dzf(z) 



6(1 -z) 



l-z 

J__z[! 

2e 2 12 



2e r[l-2e] 
ln(l - z 
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(A5) 



so we obtain the distributional identity 

(z — x)~ e 



l-z) 



-1-e 



■ 1 7T 

dx- — — = 6(1 — z) I — 

1 - x K J V 2e 2 12 



ln(l-z) 
1 - z 



(A6) 



Combining Eq.JXTJ), Eq.flMJ), and Ea.rfAH) yields the result in Eq. (EEHl) . 

To obtain the result in Eq. ()49|). we first note that we can replace g(x) with g(l) and 
then use 



fdzf{ 
Jo 



z)(l-z) 



dx 



z — x) 



(1-xf 



dz[f(z) - /(l)] / dx- ^ + f(l) / ^(1 - / dx y - 
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z - X 



dzf(z) 
plus nonsingular terms. 



1 - z 



1-xf 
(A7) 



APPENDIX B 

Here we discuss zero-bin subtractions in more detail and show how they can be imple- 
mented at the level of the Lagrangian by introducing a fictitious field. A collinear particle 
is defined in SCET by rephasing the full theory fields: 

= . (Bl) 

Loop integrals as well as phase space integrations involve both a sum over labels and an 
integral over residual momentum which can be performed using the following trick: 




which allows one to avoid doing explicit sums in loops and in phase space. Note that the 
sums in Eq. (JB1J) and Eq. (jB2J) differ. In defining the fields, we only include nonvanishing 
label momentum, p ^ 0, since when the label momentum vanishes the particle should really 
be regarded as ultrasoft. However, the trick in Eq. (jB2|) requires that the sum include the 
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so-called "zero-bin" , the momentum region with vanishing label momentum p = 0. Careful 
evaluation of the sum over labels and integrals requires modifying Eq. (jB2|) 27]: 



/ d D kF[p,k] = / d D p(F[p} - subtraction) . (B3) 

In a tree level calculation which involves phase space integrals, one can deal with the 
zero-bin by putting a cutoff in the phase space integral that separates the integration region 
into x — z ~ 0(X 2 ) and x — z ~ 0(1) regions, and only use the appropriate mode in the 
respective regions. This is somewhat clumsy and a more elegant way to acheive the desired 
result is to extend both the collinear and soft particle diagrams to cover the entire phase 
space and introduce a fictitious particle, zb , which cancels off the double counting in the 
corner of phase where such a cancellation is needed. We modify the definition of the collinear 
field so the sum is over all labels but the contribution of the spurious zero-mode is cancelled 
by a fictitious degree of freedom: 

<f>(x) = Y, e- lf *M*) ~ <A*0 • (B4) 
v 

If we use the fictitious <p zb field then we can convert sums over labels ordinary loop integrals 
and phase space intergrals using Eq. (jB2|) but there are extra diagrams with the zero-bin field 
which cancel spurious zero-bin contributions in the collinear graphs. These fields have the 
same couplings as collinear fields but their diagrams contribute with opposite sign and they 
have vanishing label momentum. Evaluating real emission graphs with (p zb yields Eq. (|6U|) . 
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